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1.  INTRODUCTION 


The  Fox- Wright  Psi  function  p'Pj.?]  and  its  normalization  are  hypergeo¬ 

metric  functions  whose  series  representations  are  given  by 
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where  ,  (z)  is  the  Gamma  function.  Thus  p'Pjz]  is  a  special  case  of  Fox’s  H-function 
H'kf  [z]  (see  e.g.  [1,  p.  50])  and  is  a  generalization  of  the  familiar  generalized 

hypergeometric  function  pFq[z]\ 
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where  the  Pochhammer  or  shifted  factorial  symbol  is  defined  by  (a)n  =  ,  (a  +  n)/,  (a) 

for  non-negative  integers  n.  Clearly,  if  we  set  A;  =  1  (i  =  1, .  . .  ,  p),  Bi  =  1 

(i  =  1, .  .  .  ,  q)  in  Equation  (1.1),  \Ii*|Y]  reduces  to  pFq[z\  given  by  Equation  (1.2). 

In  what  follows  we  shall  consider  only  the  special  case  of  the  Fox- Wright  function 
where  p  =  q  =  1;  thus 
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The  importance  of  ^^[2;]  has  recently  been  indicated  by  its  connection  with  elemen¬ 
tary  number  theory  via  Fermat’s  last  theorem  [2], [3]  and  in  applied  problems  via  the 
solution  of  trinomial  equations  [4],  Thus  it  appears  worthwhile  studying  and  record¬ 
ing  further  properties  of  the  Psi  function.  To  this  end  in  the  present  paper  we  shall 
show  that  x'PJDs;]  reduces  to  a  single  generalized  hypergeometric  function  when  A  and 
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B  assume  integer  values  and  that  when  these  same  parameters  are  rational  num¬ 
bers,  reduces  to  a  finite  sum  of  generalized  hypergeometric  functions.  Then  in 

section  4  we  shall  give  further  applications  relevant  to  the  solution  of  trinomials  in 
general  and  in  particular  in  section  6  to  certain  trinomials  which  arise  in  a  problem 
in  information  theory.  In  section  5  we  give  a  reduction  of  a  particular  case  of  Meijer’s 
G-function. 
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2.  THE  PSI  FUNCTION  FOR  INTEGER  A  AND  B 

Assume  that  A  and  B  in  Equation  (1.3)  are  positive  integers.  Clearly,  for 
A  =  B  =  1 
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so  that  in  this  case  we  have  reduction  to  the  confluent  hypergeometric  function  lEifz] 
which  converges  for  all  z  in  the  finite  complex  plane. 

By  using  Gauss’s  multiplication  theorem  for  the  Gamma  function  it  follows  that 


,  [k(n  +  //)]  =  ,  ( kfi)(kk)n(fi)n  (ji  + 

for  positive  integers  k  and  non-negative  integers  n.  Thus  since  ,  ( a  +  An)  = 
,  [A(n  +  ol/A)\  we  obtain  (cf.  [5,  p.  240,  Eq.  (1.26)]) 
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which  when  used  together  with  Equation  (1.3)  yields 
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for  positive  integers  A  and  B.  Obviously  Equation  (2.4)  reduces  to  Equation  (2.1) 
when  A  =  B  =  1. 

Similarly,  for  integer  A  >  1 
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which  follows  from  Equation  (2.3)  and  the  identity  ,  ( a  —  n)/,  (a)  =  (  —  l)n / (1  —  a)n 
with  n  replaced  by  An.  Thus  we  have  also  from  Equations  (2.3),  (2.5)  and  (1.3)  for 
A  and  B  positive  integers  the  following: 
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Note  that  the  right  hand  side  of  Equation  (2.8)  diverges  for  z  ^  0  except  when  it 
is  a  polynomial  in  which  case  it  converges  for  all  z.  In  addition,  Equations  (2.4), 
(2.6)- (2.8)  are  valid  for  either  A  =  0  or  B  =  0  by  deleting  the  parameters  which 
contain  them. 
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3.  THE  PSI  FUNCTION  FOR  RATIONAL  A  AND  B 


Assume  now  that  A  and  B  are  positive  rational  numbers.  We  may  always  assume 
that  A  =  a/ k  and  B  =  b/k  where  the  integer  pairs  a,  k  and  6,  k  need  not  be  relatively 
prime.  Thus 
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Now  partition  the  positive  integers  n  according  to  n  =  kp  —  q  where  <7  =  0, 1, .  .  . ,  &  —  1 
and  p  =  1,  2,  3, .  .  ..  Calling  for  the  moment  the  Psi  function  in  Equation  (3.1)  S  we 
obtain 
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where  the  latter  result  was  obtained  by  reversing  the  order  of  summation  by  setting 
r  =  k  —  q  in  the  penultimate  sum. 

Referring  to  Equation  (3.2)  set 
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which  we  use  together  with  Equation  (3.2)  to  obtain  for  A;  =  1,2,3,. 
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where  a  and  b  are  non-negative  integers,  the  result  being  valid  for  either  a  =  0  or 
b  =  0  by  deleting  the  parameters  which  contain  them. 

Next,  from  Equation  (2.5)  since  ,  ( a  —  An)  =  ,  \—A(n  —  a/A)],  it  easily  follows 
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where  A  is  a  positive  integer  and  n  is  a  non-negative  integer.  Equations  (2.2)  and  (3.4) 
may  then  be  used  together  with  Equation  (3.2)  to  obtain  for  non-negative  integers  a 
and  b  the  following  results  for  k  =  1,  2,  3, . .  .: 
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We  mention  here  that  obviously  results  analogous  to  Equations  (2.4),  (3.3)  etc.  may 
be  obtained  for  p'h*[2;]  defined  by  Equation  (1.1)  when  the  A;  (i  =  1, .  .  .  ,p)  and  Bi 
(i  =  1, .  .  .  ,  q)  are  integers  or  rational  numbers. 

An  interesting  corollary  may  be  obtained  from  either  Equation  (3.3)  or  (3.5)  by 
setting  a  =  6,  a  =  /3;  thus 


exp  z  = 


r  +  1 

k 


r  +  2  r  +  k 
k  k 


where  A;  is  a  positive  integer. 

We  conclude  this  section  by  noting  that  when  k  =  1,  Equations  (3.3),  (3.5)-(3.7) 
respectively  reduce  after  simplification  to  Equations  (2.4),  (2.6)-(2.8). 
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4.  APPLICATION  TO  TRINOMIAL  EQUATIONS 


Interest  in  the  solution  of  algebraic  trinomial  equations  originated  evidently  with 
Lambert  (1758).  Numerous  other  investigators  studied  them,  notably  Lagrange  (1770), 
Heymann  (1887),  Capelli  (1892),  Ramanujan  (circa  1903  and  1913),  and  Mellin 
(1915).  Berndt  [6,  p.  72,  p.  307]  gives  a  brief  history  of  the  subject  while  Belar- 
dinelli  [7,  p.  30]  presents  a  detailed  account  including  an  extensive  bibliography. 

Mellin’s  result  [8]  (see  also  [7,  p.  37]  and  [9,  p.  81])  may  be  given  elegantly  by 
employing  the  Psi  function  as  follows:  for  real  x  the  positive  root  of  the  trinomial 
equation 

yN  +  xyN~Q  -1  =  0  ,  A>Q>0  (4.1) 


is  given  by 
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Belardinelli  observed  [7,  p.  56]  that  when  A  and  Q  are  integers,  then  the  solution 
y  given  by  Equation  (4.2)  may  be  written  as  a  sum  of  A  generalized  hypergeometric 
functions  defined  by  Equation  (1.2).  However  his  result  is  incorrect.  Such  a  represen¬ 
tation  in  terms  of  generalized  hypergeometric  functions  may  be  obtained  immediately 
from  Equation  (3.7)  by  setting  a  =  1/A,  (3  =  1  +  1/A,  a  =  N  —  Q,b=Q,k  =  N, 
and  z  =  —x.  Thus  we  obtain  for  positive  integers  A  >  Q  >  1  that  the  positive  root 
of  Equation  (4.1)  is  given  by 
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where 


Z  =  {-Q)q{N  -  Q)n~Q(-x/N)n  .  (4.5) 

We  note  that  since  £  must  satisfy  |£|  <  1  in  order  for  each  jv+i-Pjv[C]  in  Equation  (4.4) 
to  converge,  then  Equation  (4.5)  yields  the  inequality  (4.3).  In  particular,  if  \x\  <  1, 
then  we  always  have  |£|  <  1  which  will  prove  useful  in  section  6. 

We  remark  that  it  is  easy  to  show  that  the  positive  solution  y  of  Equation  (4.1) 
is  also  given  by 


1 


N 


( —  + 

’  V  N  '  N1  > 


W^,(l  +  *  -  ¥r)  >'! 


■JV+I-^JV 


where 


N  N 

1  _r _ I  1  r_ _ |_  1  I  Q  —  1  r  _  1 

4  IV  T  ,v0r  JV  T  NQ  t  Q  ,  N 


N  1  NQ- 

r  +  1  r+N 

N  ’  '  '  '  ’  N 


N(N-Q)  ' 


’  N 


1  ,  N—Q  —  l 

N(N-Q)  "r  N—Q 


(4,6) 


n  =  (-Q)Q(N  -Qf-Q(x/Nf 

In  addition,  y~x  is  the  positive  root  of  zN  —  xz®  —  1  =  0. 

Other  representations  for  the  positive  root  of  Equation  (4.1)  (with  Q  replaced  by 
N  —  Q )  may  be  given.  For  by  using  [2,  Eqs.  (9)  and  (10)]  the  trinomial  equation 

yN/Q  +  xy  —  l  =  0 

for  integers  N  >  Q  >  1  has  the  positive  solution 
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thus  it  is  easy  to  see  that 
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Note  that  y~x  is  also  the  positive  root  of  zN  —  xzN~®  —  1=0. 

Further,  a  computation  similar  to  that  employed  in  deriving  Equations  (3.3)  and 
(3.5)-(3.7),  but  more  complex  in  its  details  (so  that  for  brevity  we  shall  omit  it  here), 
yields  the  result: 
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where  r]  is  given  by  Equation  (4.6),  8  is  defined  by 
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and  x  satisfies  the  inequality  (4.7). 

In  particular,  setting  Q  =  1  in  Equation  (4.8)  we  see  that  the  positive  root  of  the 
equation 

yN  +  xy  —  1  =  0,  N  >  2 


is  given  by 


jv-i  _ 

v=E(^ 


r  —  1 


—r 

~Njr-i  r! 


r  —  1 


-NFN-1 


i  _r_  _  1  v  r_  j_  N —3 

N  N- 1  ’  N">  *  *  *  ’  N  '  N- 1 

r  +  1  r-\-N  —  1 

N  ’  ‘  ‘  ‘  ’  N 


(. N  -  1  )JV“1(:r/iV)JV 
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5.  CONNECTION  WITH  MEIJER’S  G-FUNCTION 


In  [10]  Boersma  shows  for  positive  integers  a4-  (i  =  1, .  .  .  ,  p),  bi  (i  =  1, .  .  . ,  q)}  and 
k  that  the  Fox- Wright  Psi  function 


(a1,a1/k),...,(ap,ap/k)  ; 

z  =  gGi 

(Pi,b1/k),...,(Pq,bq/k)  ; 


k,a  (~Z\k  ^ 


where 


a  =  ai  ’  b  =  bt 

i= 1  i=l 

^  f\(2^brPi 

,  (Pi)---,  (Pq)  V k  M  ^ 

,  K)---,  K)  2ttV 


6  =  i-W  fe  =  i-^-tt,...,U  =  i-ai  +  ai-1 

d\  d\  d\ 

£  _1  a<2  £  _1  «2  +  l  «2  +  «2  —  1 

Sai+1  J-  5  Sai+2  J-  5  •  •  •  i  S,a\+a2 

a  2  a  2  0,2 


=  i-—  ,  £ 


+  &2  +  --  +  ap  — 1+2 


=  1  - 


+  1  +  ap  —  1 

- , ...  ,<fa  =  1 - : 


1  k  —  1 

??i  =  0  ,  ?72  =  ,rjk  = 

1  £i  _  ,  £i  +  1  1  /5i  +  —  1 

^fc+l  1  T-  ,  Vk+2  1  7  ■>•••!  kjk+bi  1  7 

»i  »i  6i 

1  £2  1  £2  +  1  1  /52  +  &2  —  1 

?7&+&i  +1  —  J-  —  T-  5  Tjk+bi  +2  —  J-  —  7  5  *  *  *  5  ?7&+&i+&2  ^  —  7 

02  62  62 


_  ,  ftq  _  fig  +  1  _  1  f3g  +  bq  —  1 

^7/c  +  6l  +&2  +  -"  +  frq-l  +1  ^  7  5  ^7&  +  &i  +62  +*"  +  ^g  —  1  +2  -L  7  5  *  *  *  5  ^//c  +  6  ^  7 

bg  bg  bg 
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In  particular,  for  p  =  q  =  1  we  get 


in 


( a,a/k )  ; 

(P,b/k)  ; 


i+b-a-fc  kb  ,  (8)  ac 

=  Z7T  2  •  ' 


a  ,  (a)  6^ 


s~ik,a 

'Ua,k+b 


V 


—z\kaa 


k  J  bb 


(5.1) 


1  a_  _  a-\-a  —  1 

a  5  *  *  *  5  a 

n  1_  k— 1  -i  j3_  -i  /3+6— 1 

"?£?•••?  A;’1  &?•••? 

Now  comparing  the  latter  result  with  Equation  (3.3)  we  obtain  for  positive  integers 


a,  6,  and  k: 


s~ik,a 

Ua,k+b 


k  )  bb 


—z  \  k  a11  |  1- 


0  i  til  1  _  £  1  _ 

'-'5  Z-.  5  *  *  *  5  k  ’  *  *  *  ’  X 


Q'-|-a  — 1 


a-b+fc-i  nr  ,  (a)  6^ 

ff+fe-i  I  V  kb  ,  (/3)  ac 

b  ' 


1  I  ’  (^)  \  "  ’  +  fcr)  fO  ci 

v  ■  (a)  k  ■  iff  +  !>')  >•!  “+1  i+* 


r  I  a  +  a  — 1 


,  /  z\k  aa 


1  r_  I  r  | 

1  ’  A:  a  ’  *  *  *  ’  k  ' 

£  I  r_  _|_  /3+6-1  r  +  l  r  +  fc  .  V&7  $ 

k  '  b  ’  *  *  *  ’  k  '  b  ’  k  ’  *  *  *  ’  A:  7 

(5.2) 


See  [11,  Ch.  V],  for  example,  for  an  introduction  to  the  G- function. 
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6.  A  PROBLEM  IN  INFORMATION  THEORY 

We  consider  a  noiseless  and  memoryless  communication  channel  [12]  with  symbols 
Si  and  s2.  The  time  for  symbol  Si  (s2)  to  pass  through  the  channel  is  the  positive 
integer  Q(N),  N  >  Q.  We  note  that  the  case  N  =  Q  is  trivial  and  will  not  be 
discussed. 

A  transmission  over  the  channel  can  be  viewed  as  a  sequence  s  whose  terms  are  Si 
or  s2.  We  define  the  length  of  s  to  be  CiQ  -\-  c2N }  where  c4-  is  the  number  of  occurrences 
of  Si  in  s.  Let  Sn,  n  an  integer,  be  the  set  consisting  of  all  sequences  of  length  n, 
and  let  |W|  denote  the  magnitude  or  cardinal  number  of  Sn.  Since  there  is  only  one 
sequence  of  length  zero  (the  empty  sequence),  |Sb|  =  1.  There  are  no  sequences  of 
negative  length,  so  |S,_|ri||  =  0,  n  ^  0. 

The  maximal  amount  of  “information,”  in  units  of  bits  per  unit  time,  that  can  be 
transmitted  over  the  channel  is  called  the  capacity  C  and  is  defined  by 

C'  =  limsup^00  log2  .  (6.1) 

The  original  definition  was  given  by  Shannon  [12,  p.  37]  who  used  the  ordinary  limit, 
which  is  not  always  defined  (i.e.,  Q  =  2,  N  =  4,  ISW+il  =  0).  The  limit  superior 
is  always  defined  since  |W|  is  bounded  from  above  by  2n.  Other  authors  have  noted 
and  corrected  Shannon’s  definition,  but  the  error  nevertheless  has  been  perpetuated 
through  much  of  the  literature.  However,  in  practice,  the  following  correct  result  [12, 
p.  37]  is  often  used  to  express  C . 

Theorem  (Shannon):  The  capacity  is  given  by  C  =  log 2y,  where  y  is  the  unique 
positive  root  of  the  equation  zN  —  zN~®  —  1=0. 

For  the  same  reason  that  Shannon’s  definition  is  flawed,  his  proof  of  the  theorem, 
using  the  asymptotic  behavior  of  finite  differences,  does  not  hold.  Certainly  others 
have  correctly  proved  the  result,  but  we  offer  below  a  novel  proof  using  complex 
analysis.  We  note  that  a  similar  method  has  been  employed  by  Kuich  [13]  in  his 
study  of  the  entropy  of  context-free  languages. 


14 


If  we  place  the  terms  si,s2,  respectively,  onto  the  end  of  nonempty  sequences  of 
length  n  —  Q}n  —  N}  we  obtain  in  each  case  a  sequence  of  length  n.  Since  the  last  term 
of  a  nonempty  sequence  of  length  n  must  be  either  si  or  s2,  we  see  that  in  general 
the  | Shi  satisfy  the  recurrence  relation 


|Sh|  —  |Sh_g|  +  |Sh_jv|  +  h( 


On 


where  S0n  =  1  if  and  only  if  n  =  0. 

OO 

We  dehne  the  ^-transform  [14]  of  a  sequence  {an}  to  be  the  power  series  'y^janzr 

n— 0 

Applying  the  ^-transform  to  the  above,  we  arrive  at  the  formal  equation 

OO  OO  OO 

E  \Sn\zn  =  ZNJ2  I  Sn\zn  +  ^E  \Sn\z11  +  1  • 


n=0 


n=0 


n=0 


The  above  result  is  valid  for  \z\  less  than  the  radius  of  convergence  of  the  power  series 

OO 

s(z)  =  E  IWE  and  so  we  have 


n=0 


s(z)  = 


1  —  z®  —  ZN 

Since  |Sh|  <  2n,  S(z)  is  analytic  in  a  neighborhood  of  z  =  0  and  so  may  be 
expressed  uniquely  as  a  Maclaurin  series.  Therefore,  the  radius  of  convergence  of 
S(z),  which  is  given  by  1/limsup^^  \j\Sn\,  is  equal  to  the  magnitude  of  the  root  of 
smallest  modulus  of  1  —  z®  —  zN  =  0. 

Now  noting  that 

1  -  \z\Q  -  \z\N  <  |1  -  -  ZN\ 

we  see  that  the  positive  root  r  of  1  —  z®  —  zN  =  0  is  the  root  of  smallest  modulus; 
see  also  [15,  p.  122,  Theorem  (27,1)].  Thus  we  have 


r  =  1/limsup^^  yjwj  ; 


and 


so 


V  =  r  1  =  limsup^^  ~S~n 
is  the  positive  root  of  zN  —  zN~®  —  1=0. 
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Next,  by  Equation  (6.1) 


C  =  limsup,^^  log2 
=  log2  limsup,^ 

=  log2  y 


and  the  theorem  is  proved. 

Thus  by  setting  x  =  —  1,  for  example,  in  Equation  (4.4)  we  have  the  following 
Corollary:  For  N  >  Q  >  1  the  capacity  C  =  C(N ,  Q)  is  given  by 

N 


where 


and 


1  JV  Ml  N~Qr) 

■2c  =  i  +  iyM''-+  r) 


N  +  -k  -  #>’)r! 


iv+l^iv 


l,/i(r),...,/r(r)  +  *L£^,v(r),...,v{r)  + 


r+ 1  r+N 

N  ’  •  •  •  ’  N 


r  1 

M0  =  T7  + 


v  [r  = 


N  N(N  -  Q)  ’  ' v  y  TV  iVQ 

£  =  (-Q)<9(iV  -  Q)n~q/Nn 


Closed  forms  for  the  capacity  for  1  <  Q  <  N  <  4  are  of  course  readily  given  since 
polynomial  equations  up  to  the  fourth  degree  are  “solvable  by  radicals.”  Thus  we 
have,  for  example,  the  following: 


2c(2,1) 


2C(3,1) 


2  0(3,2) 


2C(4’2) 


1  +  V5 


1 

3 

V3 

3 


1  +  (M)-  + 


(V 27  -  V23\1/3 
1  +  - x -  + 


5  (2  +  V20 )1/2 


VST  +  C27y2/3' 
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In  general,  however,  2C^N^  may  be  expressed  by  means  of  the  results  previously 
given  in  terms  of  a  finite  number  P  of  hypergeometric  functions  p+iFP[z]  where 
P  =  N  or  N  —  1.  In  [16]  and  [17]  extensive  and  up  to  date  collections  of  reductions 
for  single,  double,  and  multiple  generalized  hypergeometric  functions  are  provided. 
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